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^ '. Abstract 
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Recent objections (Phys.Lett. 5637, 156.) to the published Zeplin I limit (Astropart. 
Phys 23, 444) are shown to arise from misunderstandings of the calibration data 
■ and procedure, and a misreading of the data in one of the referenced papers. 
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Following our publication of a dark matter limit from the UK ZEPLIN I 
detector [1] , A Benoit et al. [2] have claimed the time constant calibration to 
be unsound and inconsistent with other published data [3] and hence that the 
limit is invalid. We consider their points in turn: 



1.1 Event separation 



Benoit et al. state their opinion that discrimination between background and 
signal can only be achieved using completely separated distributions. This ap- 
proach is fundamentally incorrect; no two distributions ever have zero overlap 
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and thus one can only ever make a statistical statement as to the degree of 
separation of the two distributions as a function of the recoil energy. Different 
techniques may need to be applied to extract the underlying nuclear recoil 
population depending on this statistical statement, the technique used in the 
ZEPLIN I analysis is detailed in our paper [1]. 



1.2 Quenching factor 

In converting our visible energy range into recoil energy Benoit et al, choose to 
use the results of Aprile et al. [4] not available at the time of the original paper. 
However it should be noted that there are now additional measurements in 
the literature [5] at even lower energy that report higher values. Although this 
discrepancy has not been resolved, the effect of varying the quenching factor 
on the dark matter limit is fully detailed in the systematics review within our 
paper. 



1.3 The time constant data of Akimov et al.[3] 

In a footnote in §6 of our paper, we draw attention to the fact that there 
is no inconsistency between [3] and our data because [3] does not provide 
nuclear recoil data below 10 keV. Benoit et al. reproduce Figure 4 of [3] but 
unfortunately do not reproduce its original caption, which states clearly that 
the points shown at 7 keV and 15 keV do not refer to those energies, but apply 
to the whole energy range from 6 to 30 keV. The Akimov et al. caption also 
explains that this wide horizontal error bar was omitted from their plot for 
clarity. For this reason those results are, as we stated in [1], are of very little 
relevance to the ZEPLIN I analysis. 



1.4 The ambient run 



Our data are reproduced by Beniot et al. in their figure 3. Our paper clearly 
states that this was shown as corroborative data, but not used in the analysis, 
since the neutron-like pulses, though identical to neutron pulses in shape, 
cannot be proved to be unambiguously from neutrons. 
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1.5 Event rates for AraBe and ambient runs 



These runs are not comparable as "source present" and "source removed". 
The AmBe runs had been intended as a preliminary calibration only, and had 
a running time two orders of magnitude shorter than that of the ambient run. 
Thus for this shorter running time the number of neutron events expected 
from the ambient population is entirely negligible; the neutrons are thus from 
the AmBe source. 



1.6 AmBe run significance 



Benoit et al. claim that visual examination of the tail of our AmBe exposure 
distribution shows no neutron events. However, this must be compared with 
the shape of the 7-ray only time constant distribution also shown in our figure. 
This is found from high-statistics 60 Co and 137 Cs Compton scattering calibra- 
tion runs to have a very sharp fall-off of the left hand branch. As discussed 
in our paper, this has been determined with high accuracy and fits a 7-ray 
distribution at all energies. We show one example of this in Figure 2 of our 
paper. A statistical analysis of the AmBe source population, together with the 
high accuracy Compton scattering runs, then shows a clearly significant ad- 
ditional population on the tail when the AmBe source is in place. From these 
data the time constant ratio can be determined with the error bar which we 
quote. Plotting the data in logarithmic or linear scale makes no difference to 
the statistical significance of the short time-constant population. 



1.7 Neutron fraction 



Benoit et al. also suggest a discrepancy in the neutron fractions in our 20- 
30 keV and 3-8 keV fits. However, although perhaps not obvious from the 
original paper, the two cases shown are from different runs with different Csl 
tagging efficiencies, as it was necessary to use different DAQ gain and tagging 
parameters to obtain data sets above and below 10 keV. The absolute event 
numbers do not, of course, affect the time constant ratio, with the statistical 
influence of the small number of nuclear recoils at low energy accounted for 
as described in our paper. 
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1.8 Efficiency of noise cut 



The energy-dependent efficiency of the photomultiplier detector noise cut was 
first estimated, from simulations, to be in the region 70-90%. As stated in our 
paper, we were able to confirm this by comparing the fitted AmBe neutron 
source response in various energy intervals below lOkeV with and without 
applying the noise cut. It is these experimental efficiencies which are used in 
our Table 2. It should also be noted that, for the second population observed 
in both the AmBe neutron source runs and ambient runs, the individual pulses 
were examined and found to be of the shape characteristic of particle events, 
and not the shape observed for electronic or photomultiplier tube noise events. 

1.9 Variations in absolute time constant 

The ambient and AmBe source runs were not in fact contemporary source in 
and source out measurements as assumed by Benoit et al, but carried out at 
different times during several months of laboratory testing, and with different 
Xe fills. Since the absolute time constant has some dependence on Xe purity, 
some variation in time constant is possible from run to run, in addition to 
the monotonic decrease of absolute time constant with energy observed in all 
runs (as noted by Benoit et al. for our high and low energy AmBe source 
runs). As stated in our paper, the intention was to perform a full underground 
neutron/gamma calibration with a similar detector configuration as for the 
underground data. However, this was curtailed due to a technical failure within 
the thermal control mechanism. 

1.10 Gamma rejection factor 

Benoit et al. express scepticism at the large background rejection factor achiev- 
able with ZEPLIN I. This good background rejection results from the fact that 
any signal would lie at the tail of the 7-ray r distribution, where the latter is 
shown by Compton scattering calibration data to be falling very sharply at all 
energies, and with entirely positive curvature in the quantity d(log N)/d(r). 
The majority of the 7-ray (and (3 particle) background lies outside the signal 
region and part of the signal region is indeed extremely well separated from 
the 7-ray population. Thus analysis is performmed using Poisson statistics 
applied to a region with very few 7-ray background events (typically - 2 in 
a 2keV energy range). This is illustrated in detail by Figures 4 and 5 of our 
paper, which show how, in a given energy range, a signal limit of 10 events 
could indeed be seen as a significant slope change at the tail of a population of 
10000 7-ray events. The significance and sensitivity then continue to improve 
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with running time. Thus, this point of concern in [2] is already discussed in 
detail in §8 of our paper. 



2 Final Comments 

It remains the case that it was not the original intention to use the prelim- 
inary laboratory neutron run as a final calibration. As stated in our paper, 
the planned high-statistics neutron calibration of ZEPLIN I was precluded by 
curtailment of the underground runs due to a technical failure which would 
have necessitated removal and significant refurbishment of the detector. Re- 
sources could not be allocated to this because of the higher priority need to 
complete the new more powerful hybrid detectors, ZEPLIN II and ZEPLIN 
III. Accordingly we have set a dark matter limit using the preliminary cali- 
bration data and a conservative time constant ratio from those data, with a 
correspondingly large error bar stated in our paper. This error bar was folded 
into the subsequent analysis to give a limit which was poorer than would have 
been the case with a more accurate neutron calibration. Our paper discusses 
in detail the impact of various systematic errors (e.g. quenching factor) on this 
final limit. More recent measurements reported by the Columbia and Coimbra 
groups [4,5] show that the quenching factor does sustain values in the region 
0.2 ±0.05 down to ~10 keV, so that this assumption now appears to be valid. 

We have shown that the objections of Benoit et al. to the ZEPLIN I dark mat- 
ter limit arise from misunderstandings of the calibration data and procedure, 
and a misreading of the data in one of the referenced papers. Accordingly, 
we believe the statistical extraction of the limit to be sound, and not over 
estimated by a factor of a 1000 as stated by Benoit et a/., with a detailed 
systematics analysis given in the appendix of our paper. 
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